clc;
clear;

left = 0;
right = 1;
N = 4;
h_partition = (right-left)/N;
basis_type = 101;
Gauss_point_number = 4;

uh = poisson_solver_1D(left,right,h_partition,basis_type,Gauss_point_number);

x = linspace(left,right,N+1);
figure;
plot(x, uh, "b");
hold on;
plot(x, exact_solution(x), "r");

xlabel('X');
ylabel('uh');
legend('Numerical Solution', 'Exact Solution');
grid on;

check_FE_solution_error_1D_linear
% check_FE_solution_error_1D_quadratic